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Abstract 

Pioneer 10 and 11 were the first probes sent to study the outer planets of the Solar System 
and Pioneer 10 was the first spacecraft to leave the Solar System. Besides their already epic 
journeys, Pioneer 10 and 11 spacecraft were subjected to an unaccounted effect interpreted as 
a constant acceleration toward the Sun, the so-called Pioneer anomaly. One of the possibilities 
put forward for explaining the Pioneer anomaly is the gravitational acceleration of the Kuiper 
Belt. In this work we examine this hypothesis for various models for the Kuiper Belt mass 
distribution. We find that the gravitational effect due to the Kuiper Belt cannot account for 
the Pioneer anomaly. Furthermore, we have also studied the hypothesis that drag forces can 
explain the Pioneer Anomaly, however we conclude that the density required for producing the 
Pioneer anomaly is many orders of magnitudes greater than the ones of interplanetary and 
interstellar dust. 

Our conclusions suggest that only through a mission, the Pioneer anomaly can be confirmed 
and further investigated. If a mission with these aims is ever sent to space, it turns out, on 
account of our results, that it will be also a quite interesting probe to study the mass distribution 
of the Kuiper Belt. 



1 Introduction 



Pioneer 10 and 11 were launched in 1972 and 1973 to study the outer planets of the Solar 
System. Both probes have followed hyperbolic trajectories close to the ecliptic to opposite 
outward directions in the Solar System [1]. Due to their simple and robust design it was 
possible to determine their position with great accuracy. During the first years of its life, the 
acceleration caused by solar radiation pressure on the Pioneer 10 was the main effect [1]. At 
about 20 AU (by early 1980s) that effect became sub-dominant and it was possible to identify 
an unaccounted effect. This anomaly can be interpreted as a constant acceleration with a 
magnitude of a = (8.74 ± 1.33) x 10~ 10 ms~ 2 and is directed toward the Sun. This effect 
became known as Pioneer anomaly. For the Pioneer spacecraft, it has been observed, at least, 
until 70 AU [1]. It was also observed in Pioneer 11 [1, 2}. 

We mention that the effect of the Pioneer anomaly on comets, asteroids and the outer 
planets has been discussed in Refs. [3, 4] and found that data is so far inconclusive. 

One of the discussed possibilities for explaining this anomaly is the gravitational effect of the 
Kuiper Belt [1]. The subject has been recently discussed [5], and claimed that if one assumes 
some fairly unrealistic assumptions for the Kuiper Belt structure, namely that it starts at about 
20 AU and has mass total mass that is about ten times greater than the usually accepted one, 
than the Pioneer anomaly can be explained. However, as we shall see, our conclusions do 
not support this claim. In this work we examine this hypothesis through the study of various 
possible models for the mass distribution of the Kuiper Belt and show that in any of them one 
does not manage to fit the constant deceleration felt by the Pioneer spacecraft as well as the 
order of magnitude of the observed acceleration. 

We also investigate the amount of matter required to explain the anomalous acceleration 
via the action of drag forces and find that, likewise for the case of the Kuiper Belt, there is a 
considerable mismatch of orders of magnitude. 

Before closing this introduction it is relevant to point out that the Pioneer anomaly requires 
for sure, some further confirmation, even though, mechanical and engineering causes seem to be 
ruled out (see however Ref. [6]). This means that evidence points toward new physics causes. 
The most promising new physics explanations for the anomaly include a scalar field with a 
suitable potential [7], the running effect of gravitational coupling [8, 9], etc. A fairly complete 
list of theoretical proposals can be found in Ref. [7]. 

This papers is organized as follows: In section 2 we discuss the main features of the Kuiper 
Belt structure, its mass distribution and the resulting gravitational acceleration. In section 3 
we consider the effect of drag forces on the spacecraft and quantify the densities required to 
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account for the anomaly. Finally in section 4, we present our conclusions. 

2 Kuiper's Belt Gravitational Force and other effects 

In this work we test the gravitational effect of four different models for the Kuiper Belt mass 
distribution: two-ring, uniform disk, non-uniform disk and uniform torus models. The non- 
uniform disk was considered in the pioneering work of Boss and Peale [10]. The two-ring and 
the uniform disk models were considered in Ref. [1]. The total mass of Kuiper Belt is taken 
to 0.3 Earth masses (this is shown not to qualitatively affect our conclusions). This is the 
maximum value allowed for the Kuiper Belt dust obtained from far-infrared emission [11]. For 
an extensive discussion on the Kuiper Belt see, for instance, Ref. [12]. 

Before examining in detail the effect of the gravitational field of the Kuiper Belt one may 
wonder whether other effects such as, for instance, the migration of debris can be at the origin 
of the Pioneer anomaly. However, evidence points otherwise. Indeed, Pluto's orbit is highly 
eccentric (e = 0.25), inclined (i = 17°) and lies in a 3 : 2 resonant orbit with Neptune. Pluto 
and plutinos have a mechanism that prevents close encounters with Neptune. The relative 
motion of Neptune and Pluto ensures that when Pluto is at a crossing point of Neptune's orbit 
they are far away. The minimum distance between them is 17 AU. Pluto has a libration of the 
argument of the perihelion, this puts the perihelion at almost the maximum distance away from 
the ecliptic, which maintains the orbital stability [13, 14, 15]. Most likely, Pluto was formed 
in a heliocentric orbit in the solar nebula instead of the planetary disk. Malhotra [15] suggests 
that Pluto was formed in a heliocentric orbit, circular and with low inclination with respect to 
Neptune. Also, it is possible that Pluto and other trans-Neptunian objects have been caught 
in resonance orbits with Neptune as a result of the formation of the outer Solar System. 

In the late stages of planet formation, planets and, in particular the giant ones, were sur- 
rounded by planetesimal's debris. Subsequently, this debris was removed by the gravity of the 
planets, thereby causing the evolution on their orbits (via energy and angular momentum ex- 
change) and the formation of the Oort Cloud. The numerical simulations performed in Ref. [16] 
show that Jupiter's orbit migrated inward, in opposition to the other giant planets orbits which 
migrated outward. Thanks to Neptune's outward migration, it has been possible to capture 
Pluto and others trans-Neptunian objects. A consequence of this evolution is an increase in the 
eccentricity and, in some cases, also an increase in the inclination of the orbit of the captured 
object. A fairly general conclusion is that, migration effects imply, for objects beyond Jupiter, 
in a pressure outward the Solar System and hence the Pioneer anomaly cannot be accounted 
by effects of this nature. 
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2.1 Mass distribution and gravitational acceleration 

Let us now consider in detail the gravitational forces generated by the various models of the 
Kuiper Belt. The models were chosen to be consistent with the main observational features of 
the Kuiper Belt. The total mass of the Kuiper Belt was considered to be the same for all the 
models and taken to be about 0.3 Earth masses, a figure that is often found in the literature 



Kuiper belt objects can be divided in three classes according to their orbits: classical, 
resonant and scattered. The first model we have studied is based on the idea that all the 
Kuiper Belt mass is concentrated in the main resonant orbits. The second considered model 
assumes that all the Kuiper Belt mass is uniformly concentrated in the ecliptic disk, likewise 
most of the classical objects. The third studied model is based on the assumption that all the 
Kuiper Belt mass is non-uniformly concentrated in the ecliptic disk as suggested for the first 
time by Boss and Peale [10]. Finally, the last model we have considered explores the idea that 
all the Kuiper Belt mass is uniformly concentrated in a torus, encompassing also the objects 
outside the ecliptic. Each of these models encompass some of the features of the Kuiper Belt. 

Actuallyour approach is the opposite of the one adopted in Ref. [5], where the Kuiper Belt 
structure and mass were chosen to fit the Pioneer anomaly. As already mentioned this implies 
that the mass of Belt is about 1.97 Earth masses, far too great to account for the observations, 
as already mentioned (0.3 Earth mass is the figure of merit [11]), and that is starts at about 
20 AU, which has little support on the facts. 

In what follows, we consider in detail the gravitational forces generated by four different 
models of the Kuiper Belt. We point out that we consider here the toroidal model that has not 
been previously discussed in the literature and also consider for the first time the gravitational 
force away from the ecliptic for all models that we discuss. 

2.1.1 The two-ring models 

We consider first the two-ring model which consists of two thin rings lying on the ecliptic and 
with radius Ri = 39.4 AU (resonance 3:2) and R 2 = 47.8 AU (resonance 2:1). The radial 
acceleration yielded by this model at a distance r from the Sun is given by: 



[11]. 




(1) 



with 



9rg((pm) 



r — RiCos6cos(<f) m — 0) 
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where <p m is the solar ecliptic longitude of the element of mass distribution, M KB is the total 
mass of the Kuiper Belt, r, 9 and are the spherical coordinates of the solar ecliptic coordinate 
system of the probe. The radial acceleration produced by the two-ring model has been obtained 
numerically and shown in Figure 1. 

2.1.2 The uniform disk model 

Next we consider the uniform disk model, which consists of an uniform hollow thin disk lying 
on the ecliptic, within distances R m in = 30 AU and R max = 55 AU. The radial acceleration at 
a distance r from the Sun is given by: 

GMkb f 27T [ Rmax 
a ud = ~-^ — / / 9ud(r m ,(f) m )dr m d(f) m (3) 



E>2 E>2 / /„ 

with 

r m (r - r m cos9cos(<p m - <P)) ( ,s 
(H + r l m - 2rr m cos9cos((f) m - <f)))* 
where r m is the solar ecliptic radius of the element of mass distribution. The radial acceleration 
produced by the uniform disk model has been obtained numerically and exhibited in Figure 2. 

2.1.3 The non-uniform disk model 

We consider now the non-uniform thin disk model discussed for the first time by Boss and Peale 
[10] and use R m in = 30 AU and Rmax = 100 AU. We have generalized the discussion of these 
authors to any point in three-dimensional space, introducing the term cos9 in their formula 
[10]. 

For this model of the Kuiper Belt the radial acceleration at a distance r from the Sun is 
given by: 



0.1GM K b f R ™* f 27T 

Q"nud 'j 

with 



D2 C>2 I r / 

x max ^min H-min •'U 



- / gnud(r m ,(f) m )d(f) m dr m (5) 



/ / \ p / v r m r - r m cos9cos((f) m - <p) 

9nud{r rn , (p m ) - J [r m) ^ 2 + ^ _ 2r m rcosecos(<f> m - 0)) 3 / 2 1 ' 

The radial acceleration produced by the non-uniform disk model has been obtained numer- 
ically and is shown in Figure 3. Analytic techniques to study the gravitational force for this 
model have been discussed in Ref. [17]. 
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2.1.4 The toroidal mass distribution model 



Finally, we consider the torus model, i.e. a toroidal mass distribution centered over the ecliptic 
with central radius R c = 42.5 AU and thickness R t = 12.5 AU. The radial acceleration at a 
distance r from the Sun produced by the Kuiper Belt is: 

GM KB r R t r 2 ^ r 2 ^ 

a tor — — 2 7r 2fi2 J J J 9tor( r m, Pmitfirn) d(j) m df3 m dr m (7) 

where, 

-(r m R c + r 2 rn cos(3 m )cos9cos((j) - <f> m ) + r m r - r m sin(3 m sin9 

9tor{r m , Pm, <Pm) ~ — " — — T~ o , D2^ ^ ' 

(h(r m , /3 m , <j> m ) + r^ + r 2 + R 2 ) 2 

and 

h(r m , (3 m , <p m ) = 2r m cos(3 m R c - 2r m rsin(3 m sin9 - (2rr m cos[3 m + 2rR c )cos6cos(<p - <p m ) (9) 

with f3 m being the angle from the ecliptic to the element of mass distribution with origin at R c . 

The radial acceleration produced by the uniform torus model has been computed numerically 
and depicted in Figure 4. 

Through the use of these expressions for the radial accelerations, we have computed the 
gravitational acceleration due to the various models of Kuiper Belt as probed by the Pioneer 
10 spacecraft. For that we have used the Pioneer data of Ref. [18]. Our results are shown in 
Figure 5. 

For all the models that we have studied, the gravitational acceleration changes sign between 
35 and 40 AU. Until about 35 AU, the acceleration is positive for all the models, while after 
40 AU, it turns negative for all models. The closest model to yield a constant acceleration, a 
distinct observational feature of the Pioneer anomaly, is the non-uniform disk model, however 
the found magnitude is much smaller than the anomalous acceleration. The model which gives 
rise to the greatest acceleration is the two-ring model. We can see that, for all the cases, 
the maximum acceleration that one obtains is 0.064 x 10~ 10 ms~ 2 , which is 0.73% of the 
observed anomalous acceleration. Even if one considers the total mass of the Kuiper Belt as 
being 1 Earth mass (Figure 6), the maximum value for the acceleration that one obtains is 
0.212 x 10~ 10 ms~ 2 , which is 2.43% of the anomalous acceleration. Hence, our results clearly 
indicate that the gravitational acceleration due to the Kuiper Belt cannot be the cause of the 
Pioneer anomaly. 

Notice that the above expressions for the acceleration have singularities in the situation 
where the craft lies inside the mass distribution. The uniform and non-uniform disk models 
have singularities for 9 ~ 0. The torus model has singularities at 9 = 3° for r inside the torus (30 
AU < r < 55 AU), at r = 35 AU for 0° < 9 < 15.75° and at r = 50 AU for 11.25° < 9 < 13.5°. 
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3 Pioneer anomaly and drag forces 
3.1 The effect of drag forces on Pioneer 10 

Having established that the gravitational acceleration of the various models of the Kuiper Belt 
cannot account for the Pioneer anomaly, we study now to which extent drag forces might have 
affected the motion of the Pioneer spacecraft. As we shall see the inferred dust distribution 
required to explain the anomaly is many orders of magnitude greater than the interplanetary 
and interstellar dust densities as well as considerably greater than ones that can be inferred 
from the Kuiper Belt that we considered in the previous section. 

In oder to perform this analysis we examine the deceleration caused by drag forces on the 
spacecraft due to the interplanetary medium. This force can be modeled as [19] : 

= _K,p(r>.(r?A, 
m s 

where p(r) is an in situ density, presumably the one of the interplanetary medium, K s an 
0(1) coefficient associated with reflection (K s = 2), absorption (K s = 1) and transmission 
(K s = 0) of particles hitting the craft, v s is the velocity of the craft with respect to the 
medium, A s the effective cross-sectional area of the craft, and m s its mass. For the Pioneer 
spacecraft v s = 11.6 — 12.2 km/s, m s = 2Alkg (after consumption of half of the fuel) and 
A s = 5.9 m 2 . Assuming that the Pioneer anomaly can be regarded as an in situ measurement 
of the acceleration, that is a s (r) = a p , the resultant density is 2.48 x 1CT 19 gcm~ 3 . We can 
compare the density profile obtained with estimated values of the interplanetary dust density 
[20, 21] 

Pipd > lO" 24 gem' 3 , (11) 

and the interstellar dust density [20] 

Pisd < 3 x 10~ 26 gem- 3 . (12) 

Density profiles of the form 

Pn{r) = Pon(r) (—) , (13) 



r 

for p 0n and r constant, have also been suggested for comparison [19]. For n = 0, the density 
is uniform and its estimated value is [19] 

Po = Poo = P P < 3 x 10~ 19 gem' 3 , (14) 

which corresponds to the one required to explain the anomaly, p p . 
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The comparison between the density obtained assuming the Pioneer anomalous acceleration 
as an in situ measurement, the interstellar and interplanetary dust densities and the uniform 
density suggested in Ref. [19] is depicted in Figure 7. It is clear that a mismatch of several 
orders of magnitude exist between the estimated amount of interplanetary and interstellar dust 
and the one required to explain the Pioneer anomaly. 

3.2 Drag forces and the Kuiper Belt 

In the previous section we have considered several models for the Kuiper Belt in order to 
estimate the resultant gravitational force acting on a spacecraft, we consider now the density 
of each of the analyzed models in order to compare with the density necessary to explain the 
anomaly due to drag forces. The total mass of Kuiper Belt is taken to 0.3 Earth masses. 

For the two-ring model, the density was given by M KB /2ir(Ri + R 2 ). We consider now the 
following three-dimensional estimate: 

Mkb ( . 

Prg ~ 2n5rSz(R 1 + R 2 ) ' [ } 

with 8z = 16 AU and 5r = 5z/2 = 8 AU. The value of Sz is estimated considering the double 
of the average of height above the ecliptic of the Kuiper Belt Objects. 

For the uniform disk model, the considered density was given by M KB /n(Rl — R\). We 
consider instead, the density in three dimensions as 

Mkb 

Pud ~ n5z(Rl - Rj) ' ( } 

On its turn, the density for the non-uniform disk model was taken to be 0.1 M KB f(r m ) / n{R\ — 
R\). In three dimensions we consider the following density 

0AM KB f(r m ) 

Pnud n5z(R 2 2 -Rl) ' 1 } 
Finally, for the uniform torus model in three dimensions is given by: 

Mkb , 1fi , 

pt = ■ (18) 

Figure 8 exhibits the density profiles of the considered Kuiper Belt models and the density 
profile obtained assuming the Pioneer anomalous acceleration as an in situ measurement. Our 
results can be summarized as follows: the densities of the two-ring model, the uniform disk 
and the toroidal model are uniform, but are about two orders of magnitude greater than the 
one necessary to explain the anomaly, p p ; for the non-uniform disk, we find that the computed 
density is greater than p p between 35 AU to about 62 AU, and smaller elsewhere. 
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4 Conclusions 



In this work we have shown that the Pioneer anomaly cannot be explained by the gravitational 
acceleration of the various Kuiper Belt models considered (two-ring, uniform and non-uniform 
disk, and toroidal mass distribution). In none of the studied cases a constant acceleration was 
found and the order of magnitude of the obtained accelerations is at best about a few percent 
of the observed anomalous acceleration. 

We have also verified that if one assumes that the Pioneer density is due to drag forces 
acting on the spacecraft the required density is many orders of magnitude greater than the 
accepted for the interplanetary and interstellar dust densities. Furthermore, we have compared 
these densities with the ones corresponding the various Kuiper Belt mass distribution models 
and found that the the densities of the two-ring model, the uniform disk and the toroidal 
model are about two orders of magnitude greater than the one necessary to account for the 
anomaly. This means that for the two-ring model, the uniform disk and the toroidal model the 
anomalous acceleration can be explained by drag forces only if about less than one percent of 
mass distribution contributes to the drag. For the non-uniform disk, explaining the Pioneer 
anomaly via a drag force can be most likely excluded as only at 35 AU and about 62 AU, the 
computed density matches p p . 

These conclusions provide further support to the idea that only through a mission (dedicated 
or not) [22, 23, 24, 25, 26, 27] this mysterious anomalous acceleration can be confirmed and more 
thoroughly scrutinized. It is quite pleasing that this possibility is being seriously considered by 
ESA in its Cosmic Vision programme for the 2015 - 2025 period [28]. 
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Figure 1: Variation of the gravitational acceleration, due to the two-ring model of the Kuiper 
Belt in terms of the distance and the solar ecliptic latitude angle. 
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(b) The full line represents r = 35AU and the dash line r = 50AU. 



Figure 2: Variation of the gravitational acceleration, generated by the uniform disk model of 
the Kuiper Belt in terms of the distance and the solar ecliptic latitude angle. 
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(b) The full line represents r = 35AU and the dash line r = 50AU. 



Figure 3: Variation of the gravitational acceleration, produced by the non-uniform disk model 
of the Kuiper Belt in terms of the distance and the solar ecliptic latitude angle. 
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Figure 4: Variation of the gravitational acceleration, generated by the uniform torus model of 
the Kuiper Belt in terms of the distance and the solar ecliptic latitude angle. 
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Figure 5: Gravitational acceleration acting on the Pioneer 10 for all the models. The lighter 
grey represents the uniform disk, the medium grey the non-uniform disk, the dark grey the 
two-ring model and the black the torus model 
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Figure 6: Gravitational acceleration acting on the Pioneer 10 spacecraft for all the models with 
a total mass for Kuiper Belt of 1 Earth mass. The lighter grey represents the uniform disk, the 
medium grey the non-uniform disk, the dark grey the two-ring model and the black the torus 
model 
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Figure 7: Density profiles of ISD (grey dash line), IPD (black full line) and the uniform density 
obtained using the Pioneer anomalous acceleration as an in situ measurement (grey full line). 
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Figure 8: Density profiles of the Kuiper Belt models with 0.3 Earth mass and the density profile 
obtained using the Pioneer anomalous acceleration as an in situ measurement. The two-ring 
model is represented by the black full line, the uniform disk by the light grey full line, the torus 
model by the black dashed line, the non-uniform disk by the dark grey full line and the the 
density obtained from Pioneer anomalous acceleration by the light grey dashed line. 
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